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In relativistic heavy ion collisions, the pr distributions of hadrons reflect the transverse motion 
generated during the collision and therefore the collision dynamics. The moderate pr values mea- 
sured to date at y/s « 20 GeV/nucleon contain a mix of hard and soft processes. The available data 
from the SPS on n, K, p, and A are used to constrain the magnitude of collective radial expansion. 
In yl + A collisions, the hydrodynamic distributions determined from the measured collective ex- 
pansion are determined up to high pr and are consistent with the 7r° data. To cleanly observe hard 
processes, measurements at pr > 4 GeV/c will be required. Pions produced in p -|- collisions 
clearly deviate from hydrodynamic distributions at pT > 2 GeV/c, underscoring the differences in 
the physics between p + A and A + A collisions. 
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The possibility of creating a new form of matter, a 
quark-gluon plasma, in ultra-relativistic heavy-ion colli- 
sions has generated intense experimental and theoretical 
effort. Transverse momentum distributions of hadrons 
reflect the transverse motion generated during the colli- 
sion and therefore the collision dynamics. If there are suf- 
ficient secondary scatterings, local thermal equilibrium 
and collective hydrodynamical motion may occur. Mea- 
sured hadronic transverse momentum (j)t) distributions 
have been used to search for such collective behavior and 
extract the hydrodynamic parameters |^. 

It is, however, the early stage of the collision, when the 
deconfined phase should exist, that is of paramount inter- 
est. Hard processes create probes for this phase, during 
the initial collisions between the participating nucleons. 
For example, particle jets arise from the fragmentation of 
high pt quarks and giuons produced in inital hard scat- 
terings. Study of high pT hadrons can therefore yield 
information about energy loss of these quarks and giuons 
as they traverse the hot, dense nuclear medium. Theo- 
retically, hard processes may be investigated using per- 
turbative QCD-inspired parton models, tuned to repro- 
duce experimental data on nucleon-nucleon (p+p) and 
nucleon- nucleus (p-l- A) collision systems . Experimen- 
tally, the challenge is to identify a pt range accessible 
with good statistical precision, but dominated by jets. 

Transverse collective expansion of the colliding system 
has been observed in high energy nucleus-nucleus (A-l-A) 
collisions ||l|, ||]. This flow boosts the pt of the hadrons, 
flattening the spectrum and effectively increasing the pr 
range where soft processes are important. The tt" pro- 
duction in Pb -I- Pb collisions at 158 GeV/nucleon has 
been successfully reproduced both by partonic models 0] 
and by hydrodynamics (i.e. soft physics dominated) 
This suggests that the complex interplay between hard 
and soft processes may not be unambiguously resolved 
by a single hadron spectrum at px < 3-4 GeV/c. 

In this paper we present a systematic study of low pr 
hadron spectra, using all available data from y/s = 17-20 



GeV/nucleon on tt, K, p and A to constrain the magni- 
tude of collective radial expansion. The expansion pa- 
rameters so obtained are used to extrapolate the pion 
spectrum to high px for each collision system. The dis- 
tribution is then compared to the data in order to de- 
termine at what Pt hydrodynamical expansion no longer 
accounts for the particle yield. As initial conditions and 
the surface profile of the expanding collision volume are 
not well constrained by the data, we choose the simplest 
profiles allowed by the data, and constrain the parame- 
ters using multiple hadronic species at transverse mass 
rriT (= — Too < 1 GeV. We include p-|-A col- 

lisions as well, because the secondary hadronic scattering 
and resulting transverse collective flow should be small. 

Collective transverse flow boosts the velocity of the 
hadrons and introduces a hadron mass dependence to the 
slope of the tot spectra. At CERN, NA44 and NA49 
observed mass dependent inverse slopes (Teff) of the 
Tot spectra, and interpreted the trend as common pro- 
duction from an expanding hadron gas. The data were 
fitted with a collective transverse flow velocity and tem- 
perature of the gas at freeze-out, when the hadrons cease 
to undergo strong interactions. The hadron spectra are 
not, however, completely exponential in tot over a large 
range, and the deviation causes the extracted velocity 
and temperature parameters to depend upon the tot fit 
range. 

We select published effective temperatures that were 
determined over a common fit range for all the data, 
(tot— Too) < 1.2 GeV, by an exponential fit (Equation |^) 
to the measured tot distributions for each hadron spec- 
trum produced at midrapidity. 
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The experiments include NA44 |l|, 0, | H and WA97 
0, |n| at the SPS; and, at the ISR, Alper et al. |l| 
and Guettler et al. For pions, the low-pT region of 

(tot — Too) < 0.3 GeV, populated by decay of baryonic 
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resonances, was systematically excluded. The effective 
temperatures used are given in Table Q with the refer- 
ences noted accordingly. 

For all but p+p and p+Be, the tt, K, p and A spectra 
were analyzed. The mass dependence of T^f / for positive 
particles is shown in Figure |l|; it is clear that a linear 
increase with particle mass exists for heavy colliding sys- 
tems. 



Equation g does not change the resulting parameters. 
We use slopes from exponential fits to the measured m,T 
spectra, as these describe the data well over the limited 
•niT range 0.3 < mx — mQ < 1.2 GeV. The simple relation- 
ship between T^ff, Tfo, and /3t was shown to agree well 
with hydrodynamical fits to the same data [Q, justifying 
the use of this simple form. 
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FIG. 1: The effective temperatures are plotted against parti- 
cle mass for positive particles produced in symmetric (A -I- A) 
collisions (left) and asymmetric (p -I- A and A -I- B) collisions 
(right) at CERN energies. Error bars are statistical only. In 
p+Pb (solid points), the effective temperature shows a slight 
dependence on particle mass. 

The dependence of T^ff on both mass and the num- 
ber of participating nucleons in the collision indicates ra- 
dial expansion. The effective temperature includes the 
local temperature of a small piece of matter and its 
collective velocity (for motion in two dimensions, the 
classical T — mo(wx)^/2 + mo(uj,)^/2 = mo{v)'^, where 
(v) = {vx) = (vy) in circular coordinates in the plane 
transverse to the collision axis). Simple exponential fits 
treat each particle spectrum as a thermal source, so the 
collective expansion velocity cannot be extracted reliably 
from a single particle spectrum; however, by using the in- 
formation from all the particles, the expansion velocity 
can be inferred. For example, NA44 ^ used 
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to separate the effective temperature into the two respec- 
tive contributions (thermal and collective motion). 

We fit equation || to the effective temperatures for par- 
ticles produced in each collision system to determine the 
conditions at freeze-out. Equation ^ assumes that the 
flow velocity (/3t) is independent of particle mass toq (all 
particles collectively travel at the same radial velocity) 
[Q. It also assumes that the pions, kaons, and protons 
have the same freeze-out temperature Tfo- 

Positive and negative hadrons yield equivalent freeze- 
out parameters. Omitting A and A from the fits with 
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FIG. 2: The freeze-out temperature Tfo (top) and mean ra- 
dial expansion velocity (/3t) (bottom) plotted for each system 
in terms of system size AB at CERN energies. The open sym- 
bols are p -I- A systems which follow the trend of increasing 
collective flow with system size. 

Figure |^ shows the dependence of the parameters on 
the size of the colliding system. The parameter Tfo re- 
mains approximately 140 MeV, regardless of system size; 
however, (Pt) increases, as would be expected from the 
larger pressure generated by copious particle production 
in the heavier colliding systems. Our values of Tfo and 
(/3t) are in reasonable agreement with other extractions 
of these parameters using two-particle correlations as well 
as single particle spectra. It should be noted that Figure 
H includes p-|-A results as well as p-|-p and A-|-A coU- 
sions (symmetric and asymmetric). The fact that p-|-Pb 
follows the trend of higher (/3t) for larger systems is re- 
markable, as significant collective flow in such collisions 
would not be expected. The effective temperatures most 
likely reflect the presence of initial state multiple par- 
ton scattering, often referred to as Cronin enhancement 
]l5[ already at relatively low hadron transverse momenta 
{pt ~ 1.3 GeV/c for pions and 1.9 GeV/c for protons). 

In heavy ion collisions, initial state multiple parton 
scattering should also be present At the relatively 
low {rriT — itiq) < 1.2 GeV fitted here, we observe a much 
larger mass dependence of the slopes in A-l-A than in 
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TABLE L Inverse slope parameters (in MeV) of hadrons for p+p, p+nucleus, S+S, S+Pb, and Pb+Pb colliding systems at 
CERN energies. The errors are statistical and systematic, respectively. 
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p+A collisions. In this soft momentum range, one would 
expect the boost from collective expansion to dominate 
over multiple scattering at the parton level. At higher 
rriT, the relative contributions should be reversed. This 
can be checked by attributing all the velocity boost to 
collective flow, extrapolating the spectrum to higher mo- 
mentum under this assumption, and comparing to data 
at higher pT or m^. 

In order to do this extrapolation, we assume that the 
soft physics processes in heavy ion collisions can be de- 
scribed by hydrodynamics. The parameters in Figure 2 
are used to extrapolate to higher transverse momentum 
the spectrum of pions from a hydrodynamically expand- 
ing gas. We use a hydrodynamics parameterization from 
U. Heinz and coworkers, which assumes that all the par- 
ticles decouple at a common freeze-out temperature Tjo, 
and that the source has an infinite length in the longitu- 
dinal direction (boost-invariance) ||l7|] . 

The transverse velocity profile of the source is parame- 
terized as a function of the radius C = i where R is the 
maximum radius of the expanding source at freeze-out, 

m)^Pse- (3) 

Ps is the maximum surface velocity, and n (= 1 for linear) 
describes the profile. A flat particle density distribution 
f{£) = const is also assumed. The geometrical average 
of the transverse velocity is therefore given by 

(ft, . Mg^i = ,4) 

Different velocity profiles, such as n = 1/2 or 2 result in 
lower or higher surface velocities, respectively. If a dif- 
ferent particle density distribution (for instance, a Gaus- 
sian) is used, then the average should be determined after 
weighting accordingly Q . Using a linear velocity profile, 
the geometrical average flow velocity in the transverse 
plane is (pt) = 2/3* /3. 



Each locally thermalized fluid element is given a trans- 
verse kick p that depends on the radial position 

p = tanh-i (A iO) ■ (5) 

The transverse mass spectrum , of the particles is 

A|rnJ(Oi^,(!^^^)/„(^^)^^ (6) 

which we integrate numerically. 

We use Equation ^ with the freeze-out temperature 
Tfo and the maximum surface velocity f3t from Figure ^ 
to extrapolate the particle spectra to high pT- The nor- 
malization A is allowed to vary, as we compare only the 
shape of the extrapolated distribution to data. Figure ^ 
shows the resulting extrapolated spectra compared to the 
neutral pion spectra in Pb+Pb, S+Pb, S+S |2§ and 
charged pions in p+W collisions (The equivalent 

p+Pb pion data in this momentum range was not avail- 
able.) The dashed lines are error bands, corresponding to 
uncertainties in Tfo and f3t from statistical and system- 
atic errors in the effective temperatures for each collision 
system. 

For heavy ion collisions, the measured ir'^ pT spectra 
are consistent with hydrodynamic distributions, though 
the highest px point tends toward the edge of the upper 
error band. The moderate pt values measured to date 
in heavy ion collisions contain a mix of hard and soft 
processes. To cleanly observe hard processes in heavy 
ion collisions, measurements at > 4 GeV/c will be 
very valuable. 

Pions produced in p+W collisions clearly deviate from 
the hydrodynamic parameterization at pt > 2 GeV/c, 
underscoring the differences in the physics between p+A 
and A+A collisions. The long tail in the p+W spectrum 
presumably reflects initial state parton multiple scatter- 
ing and the onset of hard processes, which are masked 
by the contribution of boosted soft particles in heavy ion 
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FIG. 3: The hydrodynamically extrapolated pion spectrum compared to data to the highest available pr for multiple collision 
systems at CERN. The dashed lines indicate error bands for each boosted spectrum. The hydrodynamic distribution clearly 
deviates from the p+W data at > 2 GeV/c, while no such discrepancy is observed in the heavier colliding nuclei, illustrating 
the physics differences between p+A and A+A. 



collisions. By extrapolation, the p+A data may be used 

to estimate the contribution of multiple scattering of par- 
tons to the observed broadening of the hadron spectra in 
A+A collisions. 

Stimulating discussions with Vladislav Pantucv, Xin- 
Nian Wang, Nu Xu, Axel Drees, R. Soltz, and S.C. John- 
son are gratefully acknowledged. This work was sup- 
ported by the U.S. Department of Energy under grant 
No. DE-FG02-96ER40988. 



[1] I.G. Bearden at al. (NA44 Collaboration), Phys. Rev. 

Lett. 78, 2080 (1997). 
[2] H. Appelshauser et al (NA49 Collaboration), Eur. Phys. 

J. C2, 661 (1998). 
[3] X.N. Wang, Phys. Kept. 280, 287 (1997). 
[4] X.N Wang, Phys. Rev. Lett. 81, 2655 (1998). 
[5] M.M Aggarwal ct al (WA98 Collaboration), Phys. Rev. 

Lett. 83, 926 (1999). 
[6] H. van Hecke, H. Sorge and N. Xu, Phys. Rev. Lett. 81, 



5764 (1998). 

[7] H. Boggild et al. (NA44 Collaboration), Phys. Rev. C 59, 
328 (1999). 

[8] I.G. Bearden at al. (NA44 Collaboration), Phys. Rev. C 

57, 837 (1998). 
[9] I. Bearden et al., Phys. Lett. B 388, 431 (1996). 
[10] E. Andersen et al., Journ. Phys. G 25, 171 (1999). 
[11] WA97 Collaboration, Phys. Lett. B 433, 209 (1998). 
[12] B. Alpor ct al., Nucl. Phys. B 100, 237 (1975). 
[13] K. Guettler et al., Nucl. Phys. B 116, 77 (1976). 
[14] T. Csorgo and B. Lorstad, Phys. Rev. C, 54, 1390, 
(1996). 

[15] J.W. Cronin et al., Phys. Rev. D. 11, 3105 (1975). 

[16] E. Wang and X.-N. Wang, Phys. Rev. C 64, 034901 

(2001). 

[17] E. Schnedermann, J. Sollfrank and U. Heinz, Phys. Rev. 

C 48, 2462 (1993). 
[18] S. Esumi, S. Chapman, H. van Hecke and N. Xu, Phys. 

Rev. C 55, R2163 (1997). 
[19] M.M. Aggarwal et al (WA98 Collaboration), Phys. Rev. 

Lett. 81, 4097 (1998). 
[20] R. Albrecht et al (WA80 Collaboration), Eur. Phys. J. 

C5, 255 (1998). 



